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The single transverse spin asymmetry in D meson production at RHIC can provide 
a clean measure of the gluon Sivers distribution function. At intermediate rapidity, 
D production is largely dominated by the elementary gg cc channel, where 
there cannot be any transverse spin transfer. Therefore, any transverse single 
spin asymmetry observed for D's produced in p^p interactions can only originate 
from the Sivers effect in the gluon distribution functions. A sizeable transverse 
single spin asymmetry measured by PHENIX or STAR experiments would then 
give direct information on the size of the gluon Sivers distribution function. 



Parton distribution and fragmentation functions are phenomenological 
quantities which have to be obtained from experimental observation and 
cannot be theoretically predicted. When parton intrinsic transverse mo- 
menta are taken into account, a large number of pdf's and ff's arise: the 
main difficulties in gathering experimental information on these spin and 
fc_L dependent functions is that most often more of them contribute to the 
same physical observable, making it extremely hard to estimate each sin- 
gle one separately. The Sivers function f{x,k±) ^, which describes the 
probability density of finding unpolarized partons inside a transversely po- 
larized proton, is one of these functions. It plays a crucial role since it can 
explain single spin asymmetries in terms of parton dynamics ^. 

We suggest here an experiment to be conducted at RHIC which can 
isolate the gluon Sivers effect, making it possible to reach direct independent 
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Figure 1. The unpolarized cross section for the process pp — > DX at y/s = 200 GeV, 
as a function of Ej^ and pj' at fixed pseudo-rapidity rj = 3.8 (a), and as a function of 
Xp at fixed transverse momentum pT = 1.5 GeV/c (b). 

information on the gluon Sivers distribution function A'^ fg/pr{x,k±). A 
full description can be found in a detailed recently published paper 

Let us consider the single spin asymmetry = [da'^ — da^)/ {da^ +da^) 
for p^p — > DX processes at RHIC energy, ^/s = 200 GcV. These D mesons 
originate from c or c quarks, which at LO can be created either via a 
qq annihilation, qq cc, or via a gluon fusion process, gg — > cc. The 
elementary cross section for the fusion process is much larger than the 
qq —* cc cross section. Therefore, the gluon fusion dominates the whole 
p^p — > DX process, up to xp ~ 0.6. In Fig. 1(a) wc show the unpolarized 
cross section for the process pp DX at y/s = 200 GeV as a function 
of both the heavy meson energy Ed and its transverse momentum px, at 
fixed pseudo-rapidity rj — 3.8. In Fig. 1(b) the same total cross section 
is presented as a function of xp at fixed pT = 1.5 GeV/c. The dashed 
and dotted lines correspond to the qq cc and gg — > cc contributions 
respectively. These plots clearly show the striking dominance of the gg —^ cc 
channel over most of the E]j and xp ranges covered by RHIC kinematics. 

As the gluons cannot carry any transverse spin, the elementary process 
gg cc results in unpolarized final quarks. In the qq cc process one 
of the initial partons can be polarized; however, at lowest order there is 
no single spin transfer in this s-channel interaction so that the final c and 
c are again not polarized. Consequently, the charmed quarks fragmenting 
into the observed D mesons cannot be polarized, and there cannot be any 
Collins fragmentation effect. 

Therefore, transverse single spin asymmetries in p^p DX can only 
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be generated by the Sivers mechanism, namely a spin-fe^ asymmetry in 
the distribution of the unpolarized quarks and gluons inside the polarized 
proton, coupled to the unpolarized interaction process qq cc and domi- 
nantly gg — > cc, and the unpolarized fragmentation function of either the 
c or the c quark into the final observed D meson. Detailed formulae for 
the An generated by the Sivers effect can be found in a recent paper by U. 
D'Alesio and F. Murgia ^. 

The Sivers distribution functions ^ for quarks and gluons are defined by 

A^/a/pT [Xa, k±a) = fa/pi {Xa, fe_La) - /a/pT {Xa, -k±a) , (1) 

where a can either be a light quark or a gluon. 

Here, we consider intrinsic transverse motions in the distributions of 
initial light quarks, in the elementary process and in the heavy quark frag- 
mentation function, i.e. we consider a fully non planar configuration for 
the partonic scattering. This has two main consequences: on one side, tak- 
ing into account three intrinsic transverse momenta makes the kinematics 
highly non trivial; on the other side it generates a large number of con- 
tributions, other than the Sivers effect, which originate from all possible 
combinations of k± dependent distribution and fragmentation functions, 
each weighted by a phase factor. This is a crucial topic, extensively dis- 
cussed by U. D'Alesio in these proceedings and in a paper dedicated to 
the suppression of Collins and other effects in proton-proton scattering ^. 
We have explicitely verified that all contributions to the p^p DX single 
spin asymmetry from k± dependent pdf's and ff's, aside from those of the 
Sivers functions, are multiplied by phase factors which make the integrals 
over the transverse momenta either negligibly small or identically zero. 

In the previous literature, analyses of single spin asymmetry data were 
performed under the assumption fg/pt = 0. RHIC data on in 
p^p — > DX will enable us to test the validity of this assumption. In fact, 
as the gg cc elementary scattering largely dominates the process up to 
Xp — 0.6 (sec Fig. 1), any sizeable single spin asymmetry measured in 
p^p DX at moderate xp's would be the direct consequence of a non zero 
contribution of fg/pt. 

Since we have no information about the gluon Sivers function, wc con- 
sider two opposite extreme scenarios: the first being the case in which 
the gluon Sivers function is set to zero, A^/g/pt (xa, fc_Lo) = 0, and the 
quark Sivers function A^ f^jp^ {xa, kj_a) is taken to be at its maximum al- 
lowed value at any Xa] the second given by the opposite situation, where 
A^ fg/pt = and A^ fg/pt is maximized in Xa- 
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Figure 2. Maximized values of |Ajv| for the process p^p DX as a function of E^) and 
Pt at fixed pseudo-rapidity (a), and as a function of xp at fixed transverse momentum 
(b). The dashed hne corresponds to a maximized quark Sivers function (with the gluon 
Sivers function set to zero), while the dotted line corresponds to a maximized gluon 
Sivers function (with the quark Sivers function set to zero). 



Fig. 2 shows our estimates for the maximum value of the single spin 
asymmetry in p^p — > DX. The dashed line shows I^atI when the quark 
Sivers function is set to its maximum, i.e. f^/pT (x) — 2fq/p{x), while 
setting the gluon Sivers function to zero. Clearly, the quark contribution 
to An is very small over most of the kinematic region, in both cases (a) 
and (b). The dotted line corresponds to the SSA one finds in the oppo- 
site situation, when fg/pf {x) — 2fg/p{x) and A^/^/pT = 0: here the 
asymmetry presents a sizeable maximum in the central Ed and positive x f 
energy region. This particular shape is given by the azimuthal dependence 
of the numerator of Afj ^ . 

By looking at Fig. 2 it is natural to conclude that any sizeable single 
transverse spin asymmetry measured by STAR or PHENIX experiments at 
RHIC in the region En < 60 GeV or —0.2 < xp < 0.6, would give direct 
information on the size and importance of the gluon Sivers function. 
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